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Abstract: Supersymmetric models with broken R-parity induced by lepton number vi¬ 
olating terms provide a calculable framework for neutrino masses and mixings. Within 
models with bilinear R-parity breaking six new physical phases appear which are poten¬ 
tial sources of novel CP-violating phenomena compared to the minimal supersymmetric 
extension of the standard model. We consider CP-violating observables in the decays of 
the lightest supersymmetric particle in this class of models. We show that: (i) Neutrino 
physics requires a strong correlation between three different pairs of phases, thus reducing 
the effective number of new phases to three, (ii) CP-violating phenomena in decays of the 
lightest supersymmetric particle due to new R-parity breaking phases turn out to be small, 
once constraints from neutrino physics are taken into account. We demonstrate that this 
feature does not depend on the nature of the lightest supersymmetric particle. 
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1. Introduction 

The standard model (SM) contains a source for CP violation, the complex phase of the 
Kobayashi-Maskawa matrix [Jlj]. However, it has often been argued that the SM is not able 
to explain the observed baryon/antibaryon asymmetry of the universe (see for example g) 
and thus new physics might be expected to show up in CP violating phenomena. 

In supersymmetric (SUSY) extensions of the SM, with their much larger particle con¬ 
tent, a considerable number of parameters could in principle be complex and consequently 
CP violation phenomena differ considerably from SM expectations g. It is therefore not 
surprising that the study of CP violation in the minimal supersymmetric extension of 
the standard model (MSSM) g] or its constrained version (sometimes called CMSSM or 
mSUGRA in the literature) has received quite some attention recently g g. 

The new phases present in SUSY models, however, can be restricted by existing upper 
limits on electric dipole moments [||. The general consensus is that (at least) one of the 
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following three conditions has to be realized: (i) The phases are severely suppressed |7j. 
(ii) Supersymmetric particles of the first two generations are rather heavy, with masses of 
the order of a few TeV at least (iii) There is a rather strong correlation between phases 
Q, leading to a cancellation of the different SUSY contribution to EDMs *. 

Far less work, however, has been devoted up to now to CP violation in supersymmetry 
with broken R-parity ($ p ) [Kj. R-parity breaking implies a violation of either lepton or 
baryon number. It is phenomenologically unacceptable that both types of terms are present 
|[TT| l . We will discuss only the lepton number violating terms in the following, because we 
focus on connections to neutrino physics. 

R-parity can be broken by bilinear and by trilinear terms. If both types were present 
up to (36 + 3) x 2 new parameters appear in the theory, all of which could be complex. 
A priori one could therefore expect that CP violating phenomena differ considerably from 
both the SM as well as the MSSM. 

Electric dipole moments do not constrain in a significant way the phases of either 
trilinear or bilinear terms individually, essentially since the leading contributions to the 
EDMs occur only at 2-loop level, if either only trilinear or only bilinear terms are present 
112]. t If both types of R-parity breaking terms were present contributions to EDMs 
appear at 1-loop level |l3| and thus the imaginary part of a certain product of bilinear and 
trilinear terms is more tightly constrained. However, even in the latter case the limit is not 
especially strong, considering the typical size of R-parity violating parameters expected 
from neutrino physics. It is thus fair to say that very little is known currently about the 
phases of R-parity breaking parameters. 

In this paper we study CP violation in SUSY with bilinear R-parity breaking, focusing 
mainly on aspects of those phases which are not present in the MSSM. The assumption of 
having only bilinear R-parity violation reduces the number of new phases to only six, as 
will be discussed in more detail later in the paper. 

Studies of bilinear R-parity breaking SUSY at this moment are mainly motivated 
by the recent discoveries in neutrino physics. Observations of atmospheric neutrinos by 
the Super-K collaboration [14] have confirmed the deficit of atmospheric muon neutrinos, 
especially at small zenith angles, and thus strongly point to non-zero neutrino masses and 
mixings. The preferred range of oscillation parameters from atmospheric neutrino data is 
currently (at 3er and 1 d.o.f.) ®IH 

0.3 < sin 2 6Atm. < 0.7, 1.2 x 10“ 3 eV 2 < A rn\ tm < 4.8 x 1(T 3 eV 2 . (1.1) 


Also the long-standing solar neutrino problem now provides strong evidence for neutrino 
flavour conversion, especially considering the recent measurement of the neutral current 
rate for solar neutrinos by the SNO collaboration [jU]]. If interpreted in terms of neutrino 

*The correlation “solution” is somewhat debated, see for example Q. 

U-loop contributions are proportional to the neutrino mass, see the first paper in [jl2[. 
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oscillations, the data indicate a large mixing angle between v e and v ]t — v T , with a strong 
preference towards the large mixing angle MSW solution (LMA). At 3er one has fl5|, 171 


0.25 < tan 2 0 Q < 0.83 


( 1 . 2 ) 


for 1 d.o.f., the best-fit-parameters being 


tan 2 9q = 0.44 A= 6.6 x 10 5 eV 


'0 


(1.3) 


This nicely confirms earlier hints found in Ref. [18|. 

Calculated at the 1-loop level |l9], |2(J, SUSY with bilinear R-parity breaking can 
explain the solar and atmospheric neutrino data |HJ]. However, [|u], 20] considered only 
real parameters and thus were not concerned about possible CP-violating effects. 

We will essentially follow [20| l, extending the calculation to the complex case. In the 
next section we will give mass matrices and the Higgs potential of the model. In Section 
3 we will discuss the constraints on the various phases of the model implied by current 
neutrino data. Section 4 will discuss possible CP-violating observables, before we close 
with a short conclusion. 


2. The Model 

2.1 Superpotential and soft SUSY breaking 

The supersymmetric Lagrangian is specified by the superpotential W 


W = e ab 


limuM + HYMdM + - »H a d H h u + e t LfHl 


( 2 . 1 ) 


where i,j = 1,2,3 are generation indices, a,b = 1,2 are SU( 2 ) indices, and £ is the 
completely antisymmetric 2x2 matrix, with £12 = 1. The symbol “hat” over each letter 
indicates a superfield, with Qi, Li, H d , and H u being SU( 2) doublets with hypercharges 
1, —1, —1, and 1, respectively, and U, D, and R being SU( 2) singlets with hypercharges 
— 3 , |, and 2, respectively. The couplings hu, ho and He are 3x3 Yukawa matrices, and 
// and €i are parameters with units of mass. 

Supersymmetry breaking is parameterized with a set of soft supersymmetry breaking 
terms, 


V soft = M^QTQ* + M^UiU* + M^DiD* + + M'^R,R 


Aj2 


v 2 j 


r ij2 r a* ] 


fb' 2 i 


+£ab 


+m 2 Hd H a d *H a d + m\ u H^H a u - [IM 3 A 3 A 3 + \M 2 A 2 A 2 + ^MiAiAr + h.c.] 
AZQtUjHl + A%Q\D 3 H a d + A%L\R 3 H a d - B^H a d H b u + B^Hl + h.c.l (2.2) 


In the following we assume that there is no intergenerational mixing in the soft terms. Let 


us first list the parameters which may be complex in the model defined by Eq. (2.1) and 
Eq. (2.2). Decomposed in modulus and phase these are given by e* = |e*|and ^ = 
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in Eq. (pTl|) , and A v = d\a,g{\A u \e lipA ^, | A c \e l<PA °, \A t \e lipA t}, A D = d\ag{\A d \e lipA d : 
\A s \e^ A ^ \A b \e iipA b}, A E = diag{|A e |e^, | A^e^, \A r \e^ A -}, Mi = (Mije^ 1 , M 2 = 
\M 2 \eH 2 , M 3 = |M 3 |e*^ 3 , B^i = \Bfi\e l ^ B and B^i = \Biei\e lipB i in Eq. ( |2.2| ). This means, 
in addition to the MSSM parameters the following six parameters are in general complex. 

ei = \ei\e i(fi H , B lM = . (2.3) 

As mentioned no restrictions on the size of the phases in Eq. exist; they can be 0(1). 

At this point it is appropriate to note that not all of the phases quoted above have 
a physical meaning. Any two of these CP-odd phases may be eliminated by employing 
the Peccei-Quinn- and the R symmetry U(1)pq and U(l)p, respectively. We will use this 
phase freedom later on to remove two unphysical phases. 


2.2 Scalar potential 

Next we consider the scalar potential and derive the tadpole equations. The scalar potential 
at tree level is 


y = E 


dW 


dzi 


+ Pd + Vsoft- 


(2.4) 


where Zi is any one of the scalar fields in the superpotential in Eq. dH, Vo are the D- 
terms, and V so f t is given in Eq. ( |2.2| ). The vacuum expectation values (vevs) of the Higgs 
(H u ,H d ) and lepton (Li) fields are complex in general. In order to determine the ground 
state of the potential Eq. (2.4), we use the linear parametrization 


H d = e 


w 






Hu = 


75^2 + ^ 


+ *¥>& 


B = e 


V2 


i>L + Vi 


+ iv\ 


(2.5) 


where v d , v u and Vi are positive, and we have set the phase of H u to zero since only relative 
phases are meaningful. The stationary conditions then read 

777 - = ( m H d + \l J ’\ 2 ' S j v d + v d D - \B^\v u cos(ip B + 0 ) - Vi\ii\\ei\ cos(^ + 9 - ip ei - rji) = 0 

77 -^ = (m 2 Hu + \l\ 2 + - v u D + Vi\Biei\ cos (<p Bi + Vi) ~ \B/j,\v d cos(ip B + 9) = 0 

77^ = v i D + iMlkibdCOS^ + 0 -<p ei - 7/1) + vi (^|ei| 2 + All^j + \Biei\v u cos(ip Bl + ??i) 
+u 2 e 2 ei cos(<^ ei + 7/1 - ip €2 - r/ 2 ) + u 3 e 3 ei cos(</? ei + r/i ~ <Pe 3 - %) = 0 
77J7 = v 2 D + \n\\e 2 \v d cos(p} f _ l + 9 -p e2 - V2) + ^2 (^|e2 1 2 + M| 2 ^ + \B 2 e 2 \v u cos(ps 2 + rj 2 ) 
+Ui£ie 2 cos(<y9 £2 +V2- Vex ~ Vi) + v 3 e 3 e 2 cos(^ e2 + r / 2 - ip e3 - %) = 0 
77^ = v 3 D + \^i\\e 3 \v d cos((p^ + 9 - p e3 - 773) + u 3 (je 3 | 2 + M| 3 ) + |-B 3 e 3 Kcos(</?b 3 + ?? 3 ) 
+uieie 3 cos (p t3 + 773 - ip ei - 7/1) + u 2 e 2 e 3 cos( 93 e3 + 773 - ip €2 - 7/2) = 0 ( 2 . 6 ) 
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where D = |(<? 2 + g l2 ){y\ + v 2 + v 2 + v 2 i 


v 2 ), and 


d(V) 

89 

d(V) 

drji 

d(V) 

dm 

d{V) 

dm 


v u v d \Bfi\sin(ip B + 9) + ViV d \ei\\ii\ sin (cp^ + 9 - ip €i - rji) = 0 

viVd\n\\ei\ sin(y 3 M + 9 - <p ei - r/i) + viu u |-Biei| sin(<^ Bl + r/i) 

+uir; 2 eie 2 sin(</? ei + r/i - </? e2 - m) + ^W 3 eR 3 sin(</? £1 + r}\ - ip e3 - m) = 0 
V2Vd\g\\m\ sin(^ + 9 - (p e2 - m) + ^2^1-8262! sm((p B2 + r/ 2 ) 

-■y 2 fie 2 ei sin(<^ ei + m - ip €2 - m) + v 2 v 3 e 2 e 3 sin(^ e2 + ?y 2 - <Pe 3 - m) = 0 
V3Vd\n\\m\ sin (ipfj, + 9 -<p €3 - rn) + ^3^|^3e 3 | sin((p B3 + m) 

-U 3 r;ie 3 ei sin(</? ei + m - tp e3 - m) ~~ mmmm sin(</? £2 + m ~ '~Pt 3 - m) = 0 

(2.7) 


where the repeated index i implies a summation over * = 1,2, 3. 

It was shown in Ref. |]2l|| that spontaneous CP violation does not occur in this model 
with Vi / 0, see also |22|]. In the case of explicit CP violation, where {y^, tps, ‘Pei, VtB;} / 
0,7r, the vevs are complex in general, as can be seen from Eq. (2.7). 


2.3 Scalar mass matrices 


With the solutions to Eq. ( |2.6D and (^7j) we can give the mass matrices of the matter 
fields. The mass matrix for the neutral boson fields is a 10 x 10 symmetric matrix. The 
corresponding Lagrangian is 

C S mass = -\S T M 2 s S, (2.8) 


where we have defined S T = [< 7 °, a)), #/]. 

decomposed into three 5x5 blocks: 


The mass matrix in Eq. (2.8) can be 


M| = 


M 2 ss 

M T SP 


M 2 sp 

M 2 pp 


This mass matrix is diagonalized through 


RM 2 s R t = diag(0, m 2 Sl ,..., m| g ), 


(2.9) 


( 2 . 10 ) 


where R is a orthogonal matrix. We give here only the scalar-pseudoscalar mixing block 
M 2 SP . The scalar-scalar block M 2 SS and the pseudoscalar-pseudoscalar block Mp P can 
be found in Ref. [^] for real parameters. They are valid also for the complex case after 
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obvious substitutions for complex parameters. 



0 

\Bg\s^ 

— M e i s pi 


_ M e 3|5p 3 



0 

-\B 1 e 1 \s (71 

-|-B 2 62 |s CT2 

-|-B 3 e 3 |s CT 3 

M 2 SP — 

Hki|s P i 


0 

€l l e 2 s Ai2 

|ei||e 3 |sA 13 


Hhl s P2 

— \B 2 e 2 \s a2 

— e l|l £ 2 |«Ai2 

0 

5A 23 


-l^ll e 3|s p3 

— I e 3 | s rT3 

— e l|l e 3 \S\ 13 

— |e2||e3 5a 23 

0 


( 2 . 11 ) 


where we have used the shorthand notation: = sin(0 + tp B ), s Pi = sin(y? M + 6 — ip ei — 

Vi), Soi = sin(77j + ip Bi ) and s Xij = sin (<p ei + Vi ~ Pej ~ Vj)- 

Eq. (HD together with Eq. (|2.11| ) shows that scalar-pseudoscalar Higgs transitions 
appear already at tree-level, contrary to the MSSM where this is possible only at one loop 
level. However, considering the typical size of the ^violating parameters, |ej///| < O(10 -3 ) 
], such a transition has to be expected to be rather small. 


2.4 Neutralino—neutrino mass matrix 

Here we focus on the 7x7 Neutralino-neutrino mass matrix. In the basis 
= (—*Ai, —i\ 2 , v e , z/ju, v T ), it reads 




{ M x o m T \ 

y m 0 J 


( 2 . 12 ) 


Here, the 3x4 sub-matrix m contains entries from the bilinear Ijtp parameters, 


m = 


l 2 gv x e.- irn 0 \ 

-\g'v 2 <r im \gv 2 e~ i ^ 0 |e 2 |e ^2 , 

\-yv 3 e~ ir * \gv z e0 |e 3 |e^3 ) 


and the MSSM neutralino mass matrix is given by, 


(2.13) 


M x° 


( |M, \e i( >’ 1 0 -\g'v d e- i0 

0 \M 2 \e^ \gv d e~ ie 

-\g'v d e~ ld \gv d e~ ie 0 

\ \sfv u - \gv u -|/x|e*^ 


The mass matrix _A4° is diagonalized by 


\g'v u ^ 
-\gvu 

0 ) 


(2.14) 


J\j r{y *M. {) J\f° 1 = diag(m Ui ,m. x q). (2-15) 

where (i = 1,..., 3) for the neutrinos, and (j = 1,..., 4) for the neutralinos. The method 
of perturbative diagonalization of presented in J23|] holds also in the general case of a 
complex symmetric mass matrix _A4°. 
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Only one neutrino acquires mass at tree-level and its mass is approximately given by 

Mig 2 + M 2 g' 


4 det(A4..o) 


|A|- 


where, 


h i = \e i \v d S^) +Vi \ ll \S^-^. 


(2.16) 


(2.17) 


As a result, for a realistic description of the neutrino spectrum one has to improve the 
calculation to 1 -loop order. 


A complete list of 1-loop contributions can be found, for example, in [ 201. For our 
purposes it is sufficient to consider only the b — b loop, which gives in a wide range of 
parameter space the most important contribution to the neutrino masses. The 1-loop 
corrections to the neutrino sector can be written as 


(A m L 


\bb 

)ij 


16-7T 2 


m b hl sin 2 9~ h e^A Jf 3 * A log 



(2.18) 


where h b = e ie \/2m b /v ^, 9~ b is the mixing angle and <p b is the phase in the sbottom sector, 
respectively. m~ h are the two sbottom masses. The eigenvalues up to one loop are obtained 
by diagonalizing diag(m„, ( , ra^o) + (A m v )^. The corresponding mixing matrix is AT 1 , and 
the complete mixing matrix which relates weak basis and eigenbasis, is 


AT = AT 1 -AT 0 . 


(2.19) 


The CP violating Dirac phase 6 entering the oscillation formula can be extracted by 
using the following relation 


|<5| = sin 


-i 


8 Qrn (A/25 A/"|g A/’ieA/’ig) 


cos @13 sin 2@13 sin 2 @i 2 sin 2@23 


( 2 . 20 ) 


where @ 13 = sin -1 (|A/ 35 1 ),@12 = tan -1 (|A/ 25 /A/ 15 1) and @23 = tan -1 (|A/36/A/37|)- Note that 
the quantity in the denominator of Eq. ( |2.20| ) is one of six equivalent representations of the 
Jarlskog invariant in the neutrino sector. 


3. Neutrino data and CP violating phases 
3.1 Analytical discussion 

Data from neutrino oscillation experiments require that the bilinear R-parity breaking 
parameters must obey certain conditions [2C ]. As briefly discussed in the introduction, 
atmospheric neutrino data requires A m 2 Atm to be of the order of A m? Atrn ~ C(10 -3 ) eV 2 , 
while the solar neutrino problem can be solved for Am 2 , ~ (few) 10 -5 eV 2 in the case of 
the LMA solution. $ 

^The LOW solution Aijiq ~ 10 - ' eV 2 and the quasi vacuum solution Ato|> ~ 10 -9 eV 2 can not be 
excluded at present, but are strongly disfavoured after the SNO neutral current measurement, see for 
example [M. 
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In bilinear R-parity breaking the neutrino spectrum is hierarchical and therefore m, 


and m v2 are approximately given by m 


"3 


\f^A 


Atm 


0.05 eV and m V2 — 


\/ Am ©- 


Since, as explained above, at tree level the model has only one non-zero mass eigenstate, 
it is straightforward to see that the ratio of m U2 /m U3 — ^ j Am|,/A m 2 Atni determines the 
relative importance of the one-loop mass with respect to the tree-level mass in our model. 

The hierarchy of m V2 lm V3 ~ 0.1 (in case of the LMA solution) implies that the terms 
in Eq. (2.17) must almost cancel. Thus, three pairs of phases have to obey the condition 


+ Vi - 0 - Vfi ± 7T < 0(10 2 ) . 


(3.1) 


This observation reduces the effective number of free phases in the bilinear model from six 
to three, as we will now demonstrate. 

The phases ifr have to obey the tadpole equations Eq. ( |2.7|) for an arbitrary set of input 
parameters. In the limit ip Ei + r]i — 9 — ± n = 0, the tadpole equations Eq. ( 0 ) reduce 

to 


v u v d \Bp\ sin(<^B + 9) = 0 , 
ViV u \Biei\sin((p Bi + Vi) = 0 , 


i = 1,2,3 


(3.2) 


where the equation for 6 is just the well-known MSSM relation. To find the correct min¬ 
imum of the potential, we get from Eq. (|h^) 6 = —tps and 7 ^ = ± 7 r — Eq. (0) 

now reads (p £i = ipBi ~ Pb + modulo 2n, which means that the number of indepen¬ 
dent phases is reduced to three. This leads immediately to the important result that the 
scalar-pseudoscalar mixing vanishes. This can be read of directly from Eq. ( |2.11 ) since 
every single term is zero in this limit. 

Atmospheric neutrino measurements provide an additional constraint on the bilinear 
model, due to the large angle in <-> v T oscillations, which has to be tan 2 6Atm — 

For non-zero phases this condition could only be met if 


2 a. „ |A 2 | 2 


~ 1 . 


ip e2 + r )2 ~ ±(<p £ 3 + m) ■ (3.3) 

We will investigate below to which extent Eq. ( |3.3| ) has to be satisfied numerically. 

3.2 Numerical discussion 

Before we describe our numerical analysis we specify the parameters. We remove two 
unphysical phases by applying the Peccei-Quinn symmetry U{1)pq and the R symmetry 
U(1)r. We choose ips = 0 and (p^ = n, and set all other phases of the MSSM part 
to zero. As input phases we take ip e . and iji randomly in the range [— 7 r, 7 t], and use 
the stationary conditions Eq. (2.7) to solve for the phases pBi- In order to reduce the 
number of parameters, the numerical calculations were performed in a constrained version 
of the MSSM. We have scanned the parameters in the following ranges: M 2 £ [0,1.2] TeV, 
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|^| G [0,2.5] TeV, mo G [0,1.2] TeV, Ao/mo and Bo/mo G [—3,3] and tan/3 G [2.5,10]. All 
randomly generated points were subsequently tested for consistency with the stationary 
conditions Eq. ( p+| ) and (2.7). 

Next we consider to which extent the correlations ( |3. 1[ ) have to be obeyed. With 
our choice of the phases these correlations read: <p ei + rji = 0 modulo 27 t. We study the 
allowed region of the two sums of phases ip t2 + r /2 and <p e3 + 773. From Fig. |l] we can see 
| <+ 2j3 + 7 / 2.31 <5x 10“ 3 for most of the sums ip E23 + 7/2,3- If <+ 2 + 7/2 — ±(<+ 3 + 7/3), then 
It 3 £2,3 + ? ?2,3| < 10~ 2 is possible. This feature can be traced to the constraint tan 2 9 Atm — 1 , 
see Eq. ([h]^). Fig. || confirms the additional phase correlation given in Eq. ([h]^), where we 
have plotted tan 2 9Atm versus (c p £2 +7/2)/ ((f £3 +7/3) for points in which ip e2 3 +7/2,3 > 3 x 10 -3 . 



+ V2 [rad] 

Figure 1: Phase correlation between ip e2 + 7/2 and ip e3 + 7/3 required by the constraints from 
atmospheric neutrino data. 



<+■ 2 + 7/2 

<+ 3 + 7/3 


Figure 2: Correlation between (ip e2 + 7/2)/(^ e3 + 7/3) and tan 2 9Atm ■ 

As has been pointed out previously, Mgp =0 in the limit <p £i + r/j = 0. Even if this 
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condition is not exactly fulfilled scalar-pseudoscalar Higgs mixing can be neglected in 
general provided that \ip £i + 77 *| < 10 -2 . In the sneutrino sector the situation may be 
sligthly different. Recall that in the MSSM the sneutrino is a complex scalar held. I/Lp 
splits iRe(i') and Qrn(t') by a small amount in each generation. Let us denote the associated 
mass eigenstates as z> 1,2 . Since u 1,2 are nearly degenerate large scalar-pseudoscalar mixing 
is possible in this sector also for \ip ei + 77 *| < 10~ 2 . To demonstrate this feature, let us 
consider the third generation ziy’ 2 . As an example we consider the coupling C(u^bb), 
where p) does not contain imaginary parts in the limit <p ti + rji = 0. The interaction is 
described by the Lagrangian 


c = v l T b [c 0 i bb p L + c; lhb p R ] b, 


(3.4) 


where C„i bb = + iRo i 6 ) an d R is defined in Eq. ( [2.10 ). CP invariance means 

that R 0 i 6 = 0. The extent to which this limit is saturated can be seen in Fig. |] where we 
plot (R 2 — R 2 i 6 )/(i??i 1 + R 2 1 6 ) as a function of ip £3 + 773 . For the light (dark) points the 
mass spliting is 10 ” 3 eV < \m^i — m^ 2 ] < 10 ~ 2 eV (0.1 eV < \m^,i — 777 .^ 2 1). As can be seen, 
scalar-pseudoscalar mixing in the sneutrino sector vanishes only for very small (p £i + r]i in 
this sector. 


as 

l 


a? 


+ 

H 

Hh 

Cl 

as 


1 

0.8 
0.6 
0.4 
0.2 
0 

10- 5 10 " 4 10 “ 3 10 -2 



^e 3 + Vs [rad] 

Figure 3: Scalar-pseudoscalar mixing in the nearly degenerate v R — v 1 system. For a discussion 
see text. 


4. Observables 

In this section we discuss possible CP-odd observables arising due to the new phases in 
bilinear R-parity violating couplings. Recall that these I^ p couplings are typically two to 
three orders of magnitude smaller than the R parity conserving couplings if neutrino data 
are to be explained by R-parity violation. This already implies that CP violating effects 
induced by R-parity violating parameters can at most be seen in LSP decays, because in 
all other cases either the R-parity violating branching ratios are very small or the loop 
contributions due to the couplings to R-parity conserving decay modes are tiny. 
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In the following we will discuss various possibilities for the LSP: sfermions, neutralino 
and chargino. Here we will focus on observables arising in two-body decays of the LSP, 
which are either rate asymmetries or helicity asymmetries. Note that the assumption of 
two-body decays is well motivated by the fact that LEP has not found any supersymmetric 
particle which gives a lower bound on SUSY particle masses of order 100 GeV. 

Let us first discuss a necessary condition whether a CP asymmetry is observable or not 
before considering the different LSP classes. The relevant quantity to decide whether an 
asymmetry is observable (at la), is given by x B)^ 1 , where Aqp is the CP asymmetry 
and B is the branching ratio of the decay considered. If we assume rather optimistically 
that order 10 6 LSPs can be produced at a future collider experiment, this requires that 
(A 2 CP x B)- 1 < 0( 10 6 ). 

We want to recapitulate here, that in order to construct a rate asymmetry one needs 
final-state interactions, otherwise partial decay rates are equal due to CPT invariance even 
if CP is violated. The rate asymmetry is then built through an interference of tree-level 
amplitudes and one-loop amplitudes where a pair of intermediate particles in the loop 
is on-shell. Subsequently we will discuss various one-loop diagrams giving rise to CP 
asymmetries. 


4.1 Sfermion LSP 


Consider first the case were a squark is the LSP because here the discussion is rather 
simple. The possible final states are either q' or qvj. The relevant diagrams contributing 
to possible CP asymmetries are shown in Fig. ||: a) self energy diagrams and b) vertex 
diagrams. It is obvious from this figure that in both cases all three couplings involved violate 
R-parity and, thus, the corresponding diagram gives only a tiny contribution leading to a 
negligible asymmetry. 



Figure 4: a) Self energy diagrams and b) vertex diagrams contributing to CP asymmetries of 
squark decays. 







Figure 5: a) and c) Self energy diagrams and b) and d) vertex diagrams contributing to CP 
asymmetries of slepton decays. 


For a charged slepton (£) LSP the possible final states are qq' and tv. We consider the 
following rate asymmetry, 


. T{i^£-v)-T(i^£+v) 

•Acp = —----• 

r(^ -»• £~v) + T{£ -»• t+v) 

In this case diagrams similar to those shown in Fig. ^ with appropriate replacements 
contribute. For the same reasoning as given above for the squark case their contribution 
is small. In addition to these diagrams one finds diagrams involving a VF-boson and the 
corresponding ones where the VF-boson is replaced by a charged Higgs boson as shown 
in Fig. ||. As can be easily shown the diagrams involving a VF-boson are suppressed due 
to the small fermion masses involved. In case of a self-energy contribution with a top 
quark and a bottom quark in the loop, Fig. ||a, the CP asymmetry is suppressed by the 

factor mt ™ e . For the charged Higgs boson the situation is a little bit more subtle, because 

m w ■ ■ . _ 

in the corresponding diagram of Fig. ^jc the contribution with top/bottom quarks can be 
enhanced for large tan (3. In this case the CP asymmetry is given by 


(4.1) 


»4.cp ~ 


3 Qrn{hi sin ^(gf/q, sin f3 + g^ht cos 0)g%} 
16vr M 2 


(4.2) 
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where g RR (< 74 ) are the IjL v and CP violating left and right (left) coupling of the £ to tb 
(v£). The relevant Lagrangian is given by 


C = £ b(gi P L + g[ l Pii) t + gj £Pl v + h.c. 


(4.3) 


where Pl.r = 1/2(1 75 ). The couplings can be found in [25]. In principle large tan/3 

may lead to a large CP asymmetry if £ = r. However, for large values of tan/3 the 1-loop 
contributions to the neutrino masses tend to be too large. For this reason we have not 
found any points with (A^p x B)~ l < O(10 6 ), satisfying at the same time the constraints 
from neutrino physics. 


Let us now turn to the case where the sneutrino is the LSP. It can be shown on general 
grounds that the mass splitting between 3ie(z/) and 3m(i/) is very small if the Majorana 
mass of the corresponding neutrino is tiny [^(J. Taking the known neutrino data into 
account this implies that in our case the CP-even/CP-odd mass splitting between Jfte(P) 
and 9m(ii) is typically of order eV and thus negligibly small. In the case of CP violation the 
mass eigenstates are a superposition of vf 1 and v\. As a possible CP sensitive observable 
we consider the following helicity asymmetry, 


•Thel 


ej m 


-m 




t L t l) + T ^\^ t R t R, 


(4.4) 


where £ = r, /i, e, and we sum over the two (nearly) degenerate sneutrino states, because 
they can not be resolved experimentally. 

The same classes of diagrams appear as in the case of the charged slepton decays 
discussed previously. The dominant diagrams are now the self-energy diagrams with a 
bottom quark in the loop, see Fig. f|. The contributions where a Z is exchanged in the 
s-channel (Fig. pb) is suppressed by a factor "(T," 6 . The contribution in Fig. ||a can be 
substantial for large values of tan/3. However, as was mentioned, this in turn tends to 
drive the 1-loop contribution to the neutrino masses to be too large. Moreover, there is 
a cancellation between the contributions of the two degenerate sneutrinos for z/^ e . This 
cancellation can be seen from Fig. [7| for v i where the dark points represent the result for 


i = 1 in Eq. (4.4), whereas for the light points the two contributions are summed up. 
For u e the result is the same. For zP we find that —► r + r~) ~ O(10 -5 — 10 -4 ). 

With the expected number of events, it is practically impossible to observe such a helicity 
asymmetry. 
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Figure 6: Dominant Feynman diagrams for the absorptive part of the amplitude v l -> tt. 

10“ 4 
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10“ 6 
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10 8 
10“ 9 
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[GeV] 

Figure 7: Ahei as a function of mp where the dark (light) points represent the result without 
(with) sum of the contribution of the two degenerate v l . 

4.2 Neutralino LSP 

If the Xi is the LSP the decay modes may be Xi {W£, Zi/,hi/} 
consider rate asymmetries such as 

_ r(y° -> w+n - r(xg -> w-e+) 

0 m r(x? - w+z~) + r(x5 -> w-t+) ' 

The relevant part of the Lagrangian reads 

£ = vO'f'Pn x> + bf (O l * w P R + o! xw Pl)x? W~ + irfO^Pi + Lx. (4.6) 

The relevant couplings in Eq. ([0]) can be found in (2^]. Typically these couplings obey: 
|OgH ~ \0 V R h \ » \0^l\OTV This implies that final state interactions between 
Wl and Zv will not contribute significantly to the CP asymmetry in Eq. (|4.5p . As the 
next possibility we consider a final-state interaction between Wl and hv, see Fig. | 8 |. The 


. As an example we 


(4.5) 


4 ) 
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w,z 


Figure 8: Dominant Feynman diagram for the absorptive part of the amplitude Xi 

xi->ez. 


I W and 


contribution where the W + is replaced by the charged Higgs H + , is suppressed due to 
the Yukawa coupling and therefore not relevant. Setting small lepton masses to zero, 
a straightforward calculation gives an asymmetry proportional to Am( O h O '^ w ) (see 
Appendix). Since only the left coupling 0^ w (and not the much larger right coupling 
0^ w ) appears in the denominator of Eq. ( |4.5| ) the resulting (A^ p x B ) _1 is always above 
10 6 . 


4.3 Chargino LSP 

If the Xi is the LSP the possible final states are Z£, Wv and hi. We consider the CP 
asymmetry: 

r(xr ztr) - r (xt - zi+) 


Acp r(xi^z£-) + T(xt 

The relevant part of the Lagrangian reads 


Zi+) 


(4.7) 


c = l(0\ xh P L + 0% h P R )xTh + vrO u L xw P L xiW, + i'fiO^PL + O^P^xiZ, + h.c. 

(4.8) 




The full form of the left and right couplings can be found in [20, [25|]. Scanning the param¬ 
eters over the parameter ranges as described in the previous section we find that the left 
and right couplings in Eq. (4.8) typically are: | O v J x '' v \, 1 0 ( T xh \. \(/ff 7, \ « |0^ /l | ~ |0^ xz |. 
Therefore, similar to the neutralino we will consider only a final-state interaction between 
Zi and hi. The dominant Feynman diagram is displayed in Fig. |], while the calculation 
is given in the Appendix. The result of a numerical scan is shown Fig. pa and b where we 
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Figure 9: 
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(-4q P x B) 1 for as a function of a) m x and b) |5|, respectively. 
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plot («4j~;p x B)~ l as a function of rn x and of the modulus of the Dirac phase |<5| defined in 
Eq. ( 2.2C| ), respectively. We show (A^p x B)~ l for l = e; for i = /q r the result is similar. 
As we can see in Fig. |9|a and b, (A^p x B)~ l is always above 10 6 . 


5. Conclusion 

Supersymmetric models with bilinear R-parity breaking contain six new phases compared 
to the MSSM. These phases are currently only constrained by neutrino data. We have 
found that neutrino physics requires that the phases in the Ifi v sector have to fulfill the 
relation ip ti ~ <£Bi, but are otherwise not necessarily small. This in turn implies that 
scalar-pseudoscalar mixing is vanishingly small even though the phases may be maximal 
CP violating. The only exception is if a CP-even and a CP-odd state are nearly degenerate, 
in which case the mixing between these two states can be large. 

We have discussed CP asymmetries in the decays of various possible LSPs and con- 
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eluded that all of them are unmeasurably small. Although we have considered only one 
class of CP violating observables, namely rate asymmetries, we conjecture that also other 
CP odd observables are to small to be measurable. 

On the other hand, the neutrino oscillation Dirac phase 5 is not necessarily small if 
bilinear I^ p parameters are complex. This implies that possibly the only way to determine 
whether the parameters are complex is to measure CP violation in neutrino oscillations 
themselves. 
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6. Appendix 

Here we outline the calculation of the rate asymmetry for neutralino and chargino decay. 
The calculation is carried out in the Feynman-’t Hooft gauge. In the limit where mi = 0 
and \0% h \ » \O l * h \ the dominant one-loop amplitude, corresponding to the Feynman 
diagram in Fig. ||, has the same generic structure for chargino and neutralino decays. For 
the decay x £~V, with a negative charged lepton in the final state, the amplitude reads: 

n O/vvi 2 

= j—^u{pi)[2pi p (CQ + Ci + C 2 )Pr - m x C 2 r YpPL\u(p x )£ p *(-pv)gigo—^- , (6.1) 

Ci being the Passarino-Veltman functions, go = 0^ h [0^ h ] and g r { = — C os6 w (~5 + 
sin 2 M[—^ ^+4)1 for chargino [neutralino] decay. The subscript V stands for the 
appropriate vector boson and x f° r the appropriate SUSY particle, respectively. The tree- 
level amplitude is given by: 

M (0) = iu(pihplff3 p L + 93Pr]u(Px) sP *(-Pv) , ( 6 - 2 ) 

where g% = [0^ yw ] and g§ = O[0^ w ] for chargino [neutralino] decay. The CP 

asymmetry is now 

_ r(x -*• V£~) - T(x -»• VI+) _ 2 
■ op - P( X - Vi~) + T(x - V £+) " Af(0)t^(0) 
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Am x {^smBo(2my + 3 m 2 ) + $ smCo(2my + 2 rriym x + m x + m 2 h {my — m 2 ))} 

2 — 3 my + rri^rriy 


1 2m,y %m(g%* ffp gf) 

(4vr) 2 v (\gL | 2 + |^| 2 ) 
where 

2 2 

m V — Tnt , o o N 

QmBo = 7 t—-—-—- S(m 2 y — mi) , 

1 / (m 2 — m? ) 2 \ „ „ 

StoC 0 = —vr— -j log( H-^—o-) ©("4 - , (6.4) 

\ m^/m y J 

where v = 2 mw/g and 0 denotes the step-function. From the last line in Eq. ([T)|) we can 
see that the CP asymmetry in the case of the neutralino decay is suppressed, since we have 
| <73 | << |< 73 *|. For the chargino decay this relation is reversed, as was mentioned above. 
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